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Abstract 

Snow cover over the Northern Hemisphere plays a crucial role in the Earth’s hydrology and surface energy balance, and modulates 
feedbacks that control variations of global climate. While many of these variations are associated with exchanges of energy and mass 
between the land surface and the atmosphere, other expected changes are likely to propagate downstream and affect oceanic processes 
in coastal zones. For example, a large component of the freshwater flux into the Arctic Ocean comes from snow melt. The timing and 
magnitude of this flux affects biological and thermodynamic processes in the Arctic Ocean, and potentially across the globe through their 
impact on North Atlantic Deep Water formation. 

Several recent global remotely sensed products provide information at unprecedented temporal, spatial, and spectral resolutions. In 
this article we review the theoretical underpinnings and characteristics of three key products. We also demonstrate the seasonal and spa- 
tial patterns of agreement and disagreement amongst them, and discuss current and future directions in their application and develop- 
ment. Though there is general agreement amongst these products, there can be disagreement over certain geographic regions and under 
conditions of ephemeral, patchy and melting snow. 

© 2011 COSPAR. Published by Elsevier Ltd. All rights reserved. 
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1. Introduction 

Snow covers a considerable portion of Northern Hemi- 
sphere lands during winter. It is the component of the cry- 
osphere with the largest seasonal variation in spatial extent. 
In fact accumulation and rapid melt are two of the most 
dramatic seasonal environmental changes of any kind on 
the Earth’s surface (Gutzler and Rosen, 1992; Robinson 
and Frei, 2000; Robinson et al., 1993). In the Southern 
Hemisphere, outside of Antarctica and its surrounding 
ice shelves and sea ice, snow is generally limited to smaller 
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regions such as the Andes, Patagonia and the southern 
Alps of New Zealand (Foster et al., 2008). On decadal time 
scales, snow variations over Northern Hemisphere lands 
have also been considerable (Barry et al., 1995; Brown, 
2000; Brown and Braaten, 1998; Derksen et al., 2004; Frei 
et al., 1999; Mote, 2006; Mote et al., 2005; Ye et al., 1998), 
with declines in spring associated with warmer conditions 
(Brown et al., 2010; Groisman et al„ 1994; IPCC, 2007; 
Leathers and Robinson, 1993). Recent reports on changes 
in the Arctic environment cite snow as one of the critical 
variables (ACIA, 2004; AMAP, 2011). The expectation 
during the 21st century is that changes will be increasingly 
dramatic (Frei and Gong, 2005; Raisanen, 2007; Ye and 
Mather, 1997) and spatially and temporally complex 
(Brown and Mote, 2009; Nolin and Daly, 2006). 
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While large scale changes in snow cover are useful as 
indicators of climatic variations, snow also affects other 
components of the Earth system at a variety of scales. By 
virtue of its radiative and thermal properties which modu- 
late transfers of energy and mass at the surface-atmosphere 
interface (Zhang, 2005), snow affects the overlying atmo- 
sphere (Barry, 2002; Barry et al., 2007; Cohen, 1994; Ellis 
and Leathers, 1999; Mote, 2008; Walsh, 1984) and thereby 
plays an important role in the complex web of feedbacks 
that control local to global climate. For example, because 
of the high albedo of snow, its presence can change the sur- 
face energy balance over land and ice and therefore affect 
climate (i.e. the snow-albedo feedback). Snow also modu- 
lates the hydrologic cycle (Dyer, 2008; Graybeal and 
Leathers, 2006; Leathers et al., 1998; Todhunter, 2001); 
influences ecosystem functioning (Jones et al., 2001); and 
is a significant resource for many mid latitude populations 
and for populations whose water is derived from moun- 
tainous and northerly cold regions (Barnett et al., 2005; 
Barry et al., 2007). Snow observations are critical for the 
validation of climate models (Foster et al., 1996; Frei 
et al., 2003, 2005; MacKay et al., 2006; Roesch et al., 1999). 

With regards to the freshwater flux to the ocean, the role 
of snow is to modulate seasonal timing, and in some cases 
the amount, of discharge into the oceans. While this can 
affect coastal systems across mid-latitudes, of particular 
relevance is the fresh water flux into the Arctic basin. 
The drainage area into the Arctic Ocean is ~1.5 times the 
surface area of the Arctic Ocean itself (Peterson et al., 
2002) and river runoff is the largest source of freshwater 
input into the Arctic basin (Arnell, 2005; Miller and Rus- 
sell, 2000). Much of Arctic precipitation is derived from 
snow fall, and much of the river runoff is derived from 
snow melt. During the past century, both high latitude pre- 
cipitation (Zhang et al., 2007) and river runoff to the Arctic 
basin have increased; both are expected to increase further 
in a warming climate (Peterson et al., 2002), although the 
rates of change and relative impacts on ocean circulation 
vary spatially (Rennermalm et al., 2007). 

The studies described above do not include all the pos- 
sible nonlinear feedbacks in which snow plays a role in 
the Arctic environment (Hinzman et al., 2005). For exam- 
ple, due to the insulating effect of snow cover, changes in 
the timing of snow onset or disappearance, or changes in 
the amount of snow, may influence the state of the under- 
lying permafrost, which has been warming for decades 
(Romanovsky et al., 2010) and which is expected to deteri- 
orate during this century (Lawrence and Slater, 2005) and 
may further increase the freshwater flux. Thawing perma- 
frost may also result in a significant release of carbon to 
the atmosphere as the result of microbial decomposition 
of currently frozen organic carbon (Schuur et al., 2008). 
According to Betts (2000) the expected expansion of the 
boreal forest may lead to both negative feedbacks (an addi- 
tional carbon sink) and positive feedbacks (an albedo 
decrease) on global climate, and the net effect will be a posi- 
tive feedback with increased warming. The feedbacks 


between snow, permafrost, and freshwater flux to the Arc- 
tic Ocean associated with these processes are poorly under- 
stood (Francis et al., 2009; Rawlins et al., 2010). 

While an increased freshwater flux to the Arctic has 
potential effects on thermodynamic and ecological pro- 
cesses in the coastal zone, perhaps most importantly such 
increases have been shown in the past to diminish or com- 
pletely halt the formation of North Atlantic Deep Water 
(NADW) (Rahmstorf, 2000). This occurs because freshwa- 
ter export to the North Atlantic Ocean, the region of 
NADW formation, decreases surface water density. Model 
simulations suggest that the magnitude of expected runoff 
changes during this century may approach critical thresh- 
olds for NADW formation (Arnell, 2005; Miller and Rus- 
sell, 2000; Peterson et al., 2002). In a recent study, NADW 
formation as well as permafrost degradation and changes 
to the tundra and boreal forest ecosystems (all of which 
can be affected by snow, and all of which can affect the 
freshwater flux to the ocean) have been listed among the 
potentially critical components of the Earth system that 
may be in danger of approaching “tipping points” (Lenton 
et al., 2008). Thus, accurate monitoring of high latitude 
snow remains an essential goal. 

Because of the large extent of terrestrial snow cover and 
the difficulties in obtaining ground measurements over cold 
regions, remote sensing represents an important tool for 
studying snow properties at regional to global scales. In 
recent years, advances in satellite capabilities, as well as 
in algorithm development, have led to improved monitor- 
ing of snow across the globe. The purpose of this article 
is to review the current generation of satellite-derived glo- 
bal snow observations that has become available during 
the first decade of the twenty first century, with emphasis 
on land surfaces of the Northern Hemisphere. Theoretical 
considerations for the remote sensing of snow, and key 
products are discussed. 

2. Theoretical background 

Due to the nature of interactions between snow cover 
and electromagnetic radiation of different frequencies, 
snow can be distinguished from other terrestrial surfaces 
using satellite observations based on a number of different 
active and passive techniques (Dozier, 1989; Nolin, 2010). 
The two types of instruments used for monitoring global 
scale snow variations rely on either ( 1 ) a combination of 
the visible and infrared, or (2) microwave, portions of the 
electromagnetic spectrum (Hall et al., 2005; Matzler, 
1994; Rango et al., 2000; Scherer et al., 2005; Schmugge 
et al., 2002). These methods are limited by a number of fac- 
tors, such as clouds, forest cover fraction, terrain heteroge- 
neity and precipitation. For example, interpretation of 
visible and infrared as well as passive microwave images 
can be difficult where complex terrain causes considerable 
spatial variation within each remotely-sensed footprint of 
snow depth, surface characteristics, and satellite viewing 
angles. Nevertheless, products based on these observations 
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have been vital for monitoring snow and for our under- 
standing of the role of snow in the Earth system. Though 
global active microwave data (e.g., QuikSCAT) can also 
be used to study snow extent and depth at relatively large 
spatial scale (Tedesco and Miller, 2007a, b), data are avail- 
able only from 1999 to 2009 (when the satellite failed well 
past its expected lifetime; see http://www.jpl.nasa.gov/ 
news/news. cfm?release=2009-175 downloaded November 
201 1 ). In contrast, passive microwave data have been avail- 
able since the late 1970s, and continue to be available. At 
regional scales, airborne data can also be collected before 
and after the snow falls to study the attenuation introduced 
by the snow pack on naturally emitted gamma radiation 
(Carroll, 1987). However, the data collected with this 
method have low temporal resolution (seasonal scale) and 
cannot be used for global scale studies. Consequently, we 
focus our analysis on snow parameters estimated by means 
of visible and infrared and passive microwave sensors. 

2.1. Visible and near-infrared 

Snow extent (i.e. presence or absence of snow, regardless 
of snow amount) is, in many circumstances, relatively 
straightforward to observe using visible observations 
because of the high albedo of snow (up to ~80% or more 
in the visible part of the electromagnetic spectrum) relative 
to most land surfaces. However, limitations exist. First, vis- 
ible imagery is limited to that portion of the surface illumi- 
nated by sunlight; thus darkness and low illumination are 
problematic. Second, clouds impede visible evaluation in 
two ways. All but the thinnest clouds reflect a significant 
portion of visible radiation, preventing any visible radiative 
information about the surface from reaching the satellite. 
And, because the albedos of clouds and snow are often sim- 
ilar, the discrimination between cloud-covered and snow 
covered surfaces can be difficult. However, near-infrared 
bands can be used to distinguish between snow and most 
clouds because the near-infrared reflectance of most clouds 
is high while the near-infrared reflectance of snow is low. 

Third, vegetation can obstruct visible and infrared infor- 
mation about snow from reaching the satellite sensor. For- 
est canopies protrude above the snow pack, lowering the 
surface albedo (Robinson and Kukla, 1985) and partially 
or completely obscuring the underlying surface, making it 
difficult to determine snow extent or amount (Chang 
et ah, 1996; Derksen, 2008; Goita et al., 2003; Klein 
et al., 1998; Nolin, 2004). 

Lastly, surface heterogeneity can play a role the inter- 
pretation of visible and infrared imagery in a number of 
ways. Of particular relevance to the monitoring of high lat- 
itude snow is the presence of numerous frozen lakes in Arc- 
tic regions, which may contribute to the overestimation of 
snow covered area from visible and infrared based imagery 
during periods when lakes remain frozen after the snow has 
melted on adjacent land surfaces (Derksen et ah, 2005a; 
Frei and Lee, 2010), at least when high resolution land sur- 
face data sets are not used to filter out the signal from lake 


surfaces. Passive microwave based estimates of SWE may 
be underestimated due to the presence of lakes (Derksen 
et ah, 2005a; Rees et ah, 2006). On the other hand, surface 
heterogeneity may assist in the interpretation of snow-cov- 
ered versus snow-free ground, and of snow-covered versus 
cloud-covered scenes, when trained analysts are mapping 
snow extent using visible imagery. 

2.2. Passive microwave 

Because snow grain dimensions can be similar to micro- 
wave wavelengths, snow is efficient at scattering the micro- 
wave radiation naturally emitted from the Earth’s surface 
(Matzler, 1994). Therefore, microwave emission from a 
snow covered surface is diminished relative to a snow-free 
surface, and the presence of snow can frequently be identi- 
fied (Chang et al., 1976; Grody, 2008; Hall et al., 2005; 
Matzler, 1994; Tait, 1998; Tedesco and Kim, 2006). Fur- 
thermore, because under ideal circumstances the amount 
of scattering is proportional to the number of snow grains, 
microwave instruments offer the possibility of estimating 
the mass per unit area of water in the snow pack, which 
is often measured as snow water equivalent (SWE). In con- 
trast to visible and infrared, passive microwave does not 
depend on the presence of sunlight and thus provides an 
alternative at high latitudes; and, passive microwave is lar- 
gely (but not completely) transmitted through non-precip- 
itating clouds, offering the potential to estimate snow cover 
under many cloudy conditions that preclude visible and 
infrared observations. In practice, research using passive 
microwave exploits the fact that microwave scattering by 
ice crystals is frequency-dependent: higher frequencies 
within the microwave portion of the spectrum are scattered 
more efficiently than lower frequencies, enabling the use of 
two or more frequency bands to estimate SWE (Chang 
et al., 1987; Derksen, 2008; Derksen et al., 2005b; Grody 
and Basist, 1996). Other methods have also been evaluated 
such as one based on the inversion of a snow emission 
model (e.g., Pulliainen and Hallikainen, 2001). Clifford 
(2010) provides a review of global estimates of snow water 
equivalent from passive microwave. 

Limitations to the monitoring of snow using passive 
microwave sensors are due to a variety of factors. One 
major limitation is the presence of liquid water in the snow 
pack, the microwave emission from which masks the snow 
signal and inhibits the ability of microwave sensors to 
detect wet snow. Also, because of the relatively weak 
microwave signal emitted by terrestrial surfaces, micro- 
wave sensor footprints are necessarily large (~25 km). 
Uncertainties in snow depth and SWE estimates are associ- 
ated with the physical structure of snow packs (ice lenses, 
grain size variations and vertical heterogeneity) which vary 
in space (Chang et al., 1976; Sturm et al., 1995) and time 
(Langham, 1981) and can alter the scattering and emission 
characteristics of the snow pack. Snow pack metamorpho- 
sis, which in the Arctic region typically results in a layer of 
depth hoar (with large crystal size) near the bottom of the 
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pack, results in more efficient microwave scattering. Thus, 
a signal change measured at the satellite sensor due to snow 
metamorphosis can mimic a signal change due to a change 
in SWE. Vegetation in and above the snow pack emits 
microwave radiation, and can confound any detection 
algorithm (Chang et al., 1996; Foster et al., 1997; Tedesco 
et al., 2005). 

Finally, as a snow pack reaches a certain critical depth 
the relationship between snow-amount and MW brightness 
temperature reverses (Derksen, 2008; Markus et al., 2006; 
Matzler, 1994; Schanda et al.. 1983). When SWE exceeds 
~ 1 50 mm emission by the snow pack of microwave band 
radiation is greater than scattering, resulting in a positive 
relationship between SWE and brightness temperature. 
This is an additional source of uncertainty in SWE retri- 
evals for deep snow packs. 

2.3. Remote sensing of snow in complex terrain 

Some of the difficulties inherent in the interpretation of 
remotely sensed images are exacerbated in regions with 
complex terrain (Dozier, 1989). Due to variability of slope, 
aspect, and land cover, the local solar illumination angle 
varies within one satellite footprint. In fact, due to co-reg- 
istration differences between an image and a digital eleva- 
tion model, illumination angles, and therefore reflectance 
characteristics, are often unknown. To address such issues, 
Painter et al. (2009) developed the MODSCAG model, 
which estimates mean grain size and fractional snow cover 
from MODIS data using linear spectral mixture analysis 
and a library of reflectance characteristics of different sur- 
face types. This model has relatively small errors, and 
could potentially be applied globally, but so far has been 
validated mostly in regions of complex terrain. 

A recent study of different algorithms for estimating 
SWE from passive microwave radiances in a basin with 
complex terrain in the Canadian Rockies finds that the tra- 
ditional algorithms which are based on brightness temper- 
ature differences in difference wavelength intervals (as 
discussed above) are less accurate than Artificial Neural 
Network (ANN) techniques which can be trained on obser- 
vations from surface stations (Tong et al., 2010a). Unfortu- 
nately, due to the limited distribution of stations for 
training the ANN in their test region, they are unable to 
accurately estimate spatial variations of SWE. Tong and 
Velicogna (2010) and Tong et al. (2010b), using surface sta- 
tion observations and MODIS imagery across the Macken- 
zie River Basin, determine that the minimum or threshold 
SWE value estimated from passive microwave observations 
that can be used to determine the presence/absence of snow 
varies from sub-basin to sub-basin, depending on topogra- 
phy and vegetation cover. Nevertheless they find useful 
information in remotely sensed SWE values for hydrologic 
monitoring. As these studies indicate, the estimation of 
snow characteristics in complex terrain from remotely 
sensed imagery is an important and cutting edge field of 
study. At this time, these techniques have not been incor- 


porated into global products, and are not addressed further 
in this paper. 

2.4. Comparison and evaluation of products 

When two products disagree, which is “correct?” Is 
either one of them “correct?” Two key impediments to a 
conclusive evaluation are that there is no perfect “ground 
truth,” and that the answer depends on spatial scale. The 
most obvious method of testing the veracity of remotely 
sensed (or other gridded) products is by comparison to sur- 
face reference observations. However, there exists consider- 
able contrast between surface, or in situ, and remote snow 
observations with regards to the snow pack properties that 
can be measured, their spatial and temporal resolutions 
and domains, and the methods employed to make measure- 
ments (Brown and Armstrong, 2008; Goodison et al., 
1981). 

Even in regions with surface observations, validation 
may be difficult because of the contrasting spatial scales 
of surface and remotely sensed observations (Brubaker 
et al., 2005; Chang et al., 2005). Brubaker et al. (2005) dis- 
cuss the difficulties in comparing point measurements to 
spatially integrated satellite retrievals, especially in areas 
of sparse station networks, which are typically at high ele- 
vations and northerly regions (exactly the areas where 
snow is most prevalent). They find that there is no single 
accepted method to perform validation of remotely sensed 
snow products. Chang et al. (2005) provide an informative 
review of how varying station densities and different satel- 
lite footprints are not equally spatially representative, and 
how the differences can complicate evaluations and com- 
parisons of different products. They employ geostatistical 
techniques, as suggested by Kelly et al. (2004), to quantita- 
tively define the spatial density of station observations 
required to provide sufficient information for validation 
studies. MODIS has been found to compare well with sta- 
tion based observations as well as with the National Oper- 
ational Hydrologic Remote Sensing Center products (Hall 
and Riggs, 2007), but Riggs et al. (2005) show that even 
between different versions of MODIS snow products, anal- 
yses at different spatial resolutions can provide conflicting 
results in some cases, due to both the resolution differences 
and the averaging method. 

Despite the inherent difficulties, comparative studies to 
date have drawn some useful conclusions (Armstrong 
and Brodzik, 2001; Basist et al., 1996; Bitner et al., 2002; 
Brown et al., 2007, 2010; Derksen et al., 2003; Drusch 
et al., 2004; Foster et al., 1997; Mialon et al., 2005; Mote 
et al., 2003; Romanov et al., 2002; Savoie et al., 2007; Tait 
and Armstrong, 1996). For example, evaluations of NOAA 
visible and infrared versus passive microwave products find 
more disagreement during fall and spring than during mid- 
winter, with particular differences under forest canopies, 
over complex terrain, in areas of persistent clouds, patchy 
snow, and wet snow (Armstrong and Brodzik, 2001; Basist 
et al., 1996). Over the Tibetan Plateau these products often 
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disagree year-round (Armstrong and Brodzik, 2001; Savoie 
et al., 2007). Several recent studies identify differences 
between remotely sensed products and surface observations 
over North America during the spring ablation season 
(Brown et al., 2007, 2010; Frei and Lee, 2010). 

3. Snow products 

A number of digital products that are based on remote 
observations are available. The two visible and infrared 
based suites of products that are most widely used for 
large-scale climate research are from: (1) the Interactive 
Multisensor Snow and Ice Mapping System (IMS) (Section 
3.1) and (2) the suite of products derived from the Moder- 
ate Resolution Imaging Spectroradiometer (MODIS) (Sec- 
tion 3.2). IMS is the most recent version of a product 
that dates back to the 1960s (Matson and Wiesnet, 1981). 
IMS mapping of snow extent has relied primarily on visible 
and near infrared imagery, but includes data and informa- 
tion from a number of sources. As discussed in more detail 
below, the key feature that distinguishes IMS from other 
products is human involvement in the analysis, which is 
required for operational purposes. 

The MODIS instrument, which is used to observe a 
number of geophysical variables including snow, flies on 
NASA’s Earth Observing System (EOS) Terra satellite, 
launched in 1999. A near-twin MODIS instrument is also 
flying on board the Aqua platform, which was launched 
in 2002 (Aqua and Terra have afternoon and morning 
equatorial crossing times, respectively). Aqua also hosted 
the Advanced Microwave Scanning Radiometer - Earth 
Observing System (AMSR-E) until its failure in October 
2011. AMSR-E measured the naturally emitted radiation 
in the microwave region at five different frequencies (6.9, 
10.7, 18.7, 23.8, 36.5 and 89 GHz) at both vertical and hor- 
izontal polarizations. 

The IMS and MODIS snow algorithms both rely pri- 
marily on near-polar orbiting satellites, from which daily 
images are available at high latitudes. Other algorithms 
that have been suggested (Romanov et al., 2003; Siljamo 
and Hyvarinen, 2011) rely on geostationary satellites, 
which have the advantage of higher temporal resolution, 
but have poor spatial resolution. 

3.1. visible and near infrared based products 

3.1.1. The Interactive Multisensor Snow and Ice Mapping 
System (IMS) 

The data set that has historically been the most widely 
used for the operational mapping and climatological anal- 
ysis of large-scale snow extent (not depth or water equiva- 
lent) was produced by the US National Oceanic and 
Atmospheric Administration (NOAA) National Environ- 
mental Satellite and Data Information Service (NESDIS), 
but has been transferred to the National Ice Center 
(NIC), which is jointly supported by NOAA, the US Navy, 
and the US Coast Guard. This product has been based pri- 


marily on visible and near infrared observations, and cov- 
ers the period from late 1998 to present (Ramsay, 1998), 
with the precursor maps beginning in 1966, constituting 
the longest remotely sensed environmental time series that 
has been derived in a near-consistent fashion (Helfrich 
et al., 2007; Matson and Wiesnet, 1981; Robinson et al., 
1993). The term near-consistent is used because, due to 
changing operational requirements and evolving technical 
capabilities, this product has undergone, and continues to 
undergo, improvements and refinements (Helfrich et al., 
2007; Ramsay, 1998, 2000) as summarized briefly here. 
The two reasons for this product’s importance - as opera- 
tional input into atmospheric forecast models and as a 
long-term climatic record - are also discussed. 

Although a number of improvements and corrections in 
the production of the NOAA product occurred in the ear- 
lier years, the biggest methodological change was imple- 
mented in the late 1990s. Until that time, NOAA snow 
maps were produced on a weekly basis by trained meteo- 
rologists who would visually interpret photographic copies 
of visible band imagery, and manually produce maps that 
would subsequently be digitized with spatial resolution 
between 150 km and 200 km. In 1997 NOAA began pro- 
ducing snow maps using the IMS, with improved spatial 
(24 km) and temporal (daily) resolutions. IMS is operated 
by trained analysts who produce a daily digital product uti- 
lizing Geographic Information System technology and 
incorporating a variety of, and an ongoing expansion of, 
technological capabilities as well as sources of information. 
Since 1999, when weekly manual mapping was discontin- 
ued, daily IMS maps have been produced (Ramsay, 1998; 
Robinson et al., 1999). Technological advancements since 
1999 have led to even higher resolution (4 km) snow map- 
ping (Helfrich et al., 2007). 

IMS produces estimates of snow extent across the globe 
every day, regardless of the presence of clouds. This is pos- 
sible primarily for two reasons. First, analysts use sources 
of information other than visible and near infrared imag- 
ery. Second, because IMS analysts can loop through 
sequential images, their ability to evaluate scenes is based 
on an integration of information from both spatial and 
temporal perspectives. Thus, a key feature of the IMS 
product is that human judgment as to which data sources 
are most reliable in different conditions and regions, and 
as to the final evaluation of where the snow is, remains 
an integral part of the process, and one of the strengths 
of the IMS product. IMS also includes sea ice extent, which 
is not discussed in this report. Fig. 1 shows an example of a 
daily IMS snow map in its original projection. 

It is difficult to optimize this product for both of its two 
main uses. Its primary purpose is to provide input to atmo- 
spheric forecast models. For this purpose, continued 
product improvements are advantageous. As a record for 
evaluating long term environmental change, however, the 
value of any product is diminished if methodological 
changes (including those that provide more accurate maps) 
result in temporal inconsistencies in the data set that might 
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be difficult to distinguish from actual variations in snow 
extent. To maintain product continuity and a viable long- 
term record, IMS continues to produce a coarse (24 km) 
resolution version of the data set. And, in collaboration 
with NOAA, the Rutgers University Global Snow Lab is 
producing a climate data record in which inconsistencies 
between the earlier maps and the IMS product (in addition 
to inconsistencies during the weekly map era) are accounted 
for, and can therefore be used for ongoing analyses of his- 
torical variations (Robinson and Estilow, 2008). 

3.1.2. The Moderate Resolution Imaging Spectroradiometer 
(MODIS) 

The MODIS sensor measures radiation in 36 spectral 
bands, including the visible, near infrared, and infrared 
parts of the electromagnetic spectrum. The suite of 
MODIS snow cover products, available since 2000, are 
derived using a fully-automated algorithm that provides 
good spatial resolution (500 m), cloud detection, and fre- 
quent coverage (daily at mid to high latitudes) (Hall and 
Riggs, 2007; Hall et al., 1995, 2002; Riggs et ah, 2006). 
The MODIS snow-mapping algorithm uses a normalized 
difference between MODIS band 4 (5.45-5.65 mm) and 6 
(1.628-1.652 pm) and many additional spectral and thresh- 
old tests. In forested areas the threshold is changed based 
on results of a canopy reflectance model, using both the 
Normalized Difference Vegetation Index (NDVI) and 


Normalized Difference Snow Index (NDSI) (Klein et ah, 
1998). A thermal mask is also included to remove errone- 
ous “snow” over locations where the presence of snow is 
considered to be implausible. See Riggs et al. (2006) and 
Riggs and Hall (in press) for a description of the algorithm 
and recent upgrades. 

NASA provides a hierarchy of snow products based on 
MODIS observations, designed to satisfy the needs of a 
variety of users (http://modis-snow-ice.gsfc.nasa.gov/). 
These include a Level-2 swath product; daily and 8-day 
composite Level-3 tile products which are mapped onto 
a sinusoidal projection and available in 10° lat/lon tiles; 
as well as daily, 8-day composite, and monthly Level-3 
products available in the Climate-Modeling Grid (a lati- 
tude-longitude grid) at 0.05° resolution (Hall et ah, 2002, 
2005; Riggs et al., 2005, 2006). An 8-day composite is con- 
sidered useful because in many regions, particularly at 
high latitudes, persistent cloudiness limits the number of 
days available for surface observations (see results sec- 
tion). The Climate-Modeling Grid was developed to be 
useful for the evaluation of climate models and for studies 
at large spatial scales. Other features of the MODIS snow 
product suite include daily snow albedo (Klein and Stro- 
eve, 2002) and fractional snow cover (Salomonson and 
Appel, 2004). In addition, a new cloud-gap filled product 
provides information on cloud persistence, and uses obser- 
vations from prior days to map snow (Hall et ah, 2010). 
Fig. 2 shows examples of MODIS snow cover maps in 
swath format (MOD10_L2) following a major snowstorm 
in the northeastern United States in December 2010. The 
analysis presented in Sections 4 And 5 of this article uses 
the MOD10C1 daily Level 3 global .05° daily Climate- 
Modeling Grid product with a spatial resolution of 
~5 km. 

Validation of the suite of MODIS snow cover products 
has been undertaken by many authors as described in Hall 
and Riggs (2007). These products have also been used 
extensively in modeling efforts, both at the regional and 
hemispheric scales. A bibliography of papers utilizing the 
MODIS snow cover products for both validation and mod- 
eling may be found at: http://modis-snow-ice.gsfc.nasa.- 
gov/?c=publications . 

3.2. Passive microwave based products 

Historical passive microwave measurements are avail- 
able from the Scanning Multichannel Microwave Radiom- 
eter (SMMR) instrument (1978-1987), and the Special 
Sensor Microwave/Imager (SSM/I) instrument (1987 
through present) although some compatibility issues 
between the two products exist, due to slight differences in 
the frequency bands measured, overpass time, swath width, 
native footprint resolution, and coverage issues related to 
SMMR being powered down every other day. (Armstrong 
and Brodzik, 2001; Brodzik et al., 2007; Derksen and 
Walker, 2003). AMSR-E, available from 2002 to October 
2011, provides a suite of measurements to make it spectrally 
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Fig. 2. Images acquired on 28 December 2010 showing extensive snow cover as the result of a major snowstorm in the northeastern United States on 26- 
27 December. The image on top is a "true-color” image from the MOD09 reflectance product; the center image is a binary MOD10_L2 swath snow map; 
and the image on the bottom represents fractional snow cover from MOD10_L2, so that various shades of blue indicate different fractions of snow cover 
within each pixel. Image processing by Jeff Miller, Wylie Systems, and NASA/GSFC. 


compatible with both SMMR and SSM/1 at higher spatial 
resolution (Derksen et al., 2005b; Kelly et al., 2004). Due 
to the inherent difficulties and regional variations in the 
interpretation of passive microwave signals, the production 
of a data set that is consistently accurate across all Northern 
Hemisphere regions requires either ( 1 ) a physical approach, 
which includes robust representations of snow pack pro- 


cesses and their parameterization in retrieval schemes 
(Pulliainen and Hallikainen, 2001), or (2) a regional 
approach, which includes regionally-tuned algorithms (Fos- 
ter et al., 1997) that statistically represent regional snow 
pack processes but which are not applicable in different 
snow accumulation regimes. The physical approach is very 
challenging yet has the potential of being widely applicable 
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as our knowledge of, and ability to model, regional snow 
pack processes improves, while the statistical approach is 
applicable only in the few regions for which retrieval 
schemes have been calibrated. 

The global AMSR-E SWE product suite (Tedesco et al., 
2011 updated) consists of daily, pentad (five-day) maxi- 
mum and monthly average SWE estimates that together 
comprise the only NASA satellite-based SWE product 
available to the scientific community. As an example. 
Fig. 3 shows the snow depth obtained from the AMSR-E 
product for January 30th, 2005. 

The AMSR-E SWE operational algorithm takes advan- 
tage of several AMSR-E channels that are unavailable 
from SSM/I and SMMR. Snow depth is derived as a com- 
bination of microwave brightness temperature differences 
at different frequencies, weighted by coefficients derived 
from the difference between vertical and horizontal polar- 
izations. These coefficients replaced a previously used static 
coefficient to attempt to capture the spatio-temporal vari- 
ability of parameters such as grain size (Kelly, 2009; Kelly 
et al., 2003; Tedesco and Narvekar, 2010). The algorithm 
uses a spatially varying but temporally static map of snow 
density. 

Environment Canada (http://ccin.ca/cms/en/socc/snow/ 
swe/currentSnow.aspx) also produces a regional passive 
microwave SWE product for central Canada, including 
the Prairies and part of the boreal forest back to 1978. 
Until December 1999, this product relied on a single algo- 
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Fig. 3. Snow depth on January 30, 2005 obtained from the AMSR-E 
product. 


rithm that was calibrated for the prairies region (Goodison 
and Walker, 1995). Since that time algorithms that correct 
for the effects of different forest types (Derksen, 2008; 
Goita et al., 2003) and the sub- Arctic tundra (Derksen 
et al., 2010) have been developed. 

3.3. Combined products 

One promising avenue, and an area where great efforts 
are currently being made, is to refine our abilities to com- 
bine ground observations and models with space borne 
remotely sensed data. Tait et al. (2000) provide a helpful 
review, and describe a prototype of a fully automated prod- 
uct that includes station observations as well as both visible 
and microwave retrievals. Elere we briefly review some 
products that include combinations of satellite, station 
observations, and models. While not exhaustive, it provides 
examples of the variety of integrative products and meth- 
ods that have been produced. 

3.3.1. The Canadian Meteorological Center snow product 
( station observations and modeling ) 

Snow depth from the Canadian Meteorological Center 
Daily Snow Depth Analysis includes a hybrid modeling/ 
observational approach based on optimal-interpolation of 
daily snow depth observations from hundreds of stations 
globally, with snow density estimated from a simple snow 
pack model (Brasnett, 1999). This model output is consid- 
ered most dependable over regions with significant station 
coverage, which is generally south of 55° North, where 
model results are well constrained by observations. Over 
most of the Arctic, in contrast, where there are few obser- 
vations, the analysis is based mostly on model results, and 
is skewed towards snow depth observations at coastal loca- 
tions with observing sites at open areas near airports. Snow 
at these sites tends to be shallower and to melt out earlier 
than snow in surrounding terrain. Nevertheless, this analy- 
sis is considered to be a reasonable estimate of snow depth 
over data-poor Arctic regions, and has been used in a num- 
ber of studies (Brown and Mote, 2009). Here we use CMC 
modeled snow depths for comparison with AMSR-E snow 
depths. 

3.3.2. GlobSnow (satellite, station, and model) 

In 2008 the European Space Agency embarked on an 
effort to develop a long term snow cover data set called 
GlobSnow with sufficient homogeneity to be acceptable 
for climate change analysis. The GlobSnow product cur- 
rently includes global gridded information on snow extent 
and SWE across the Northern Hemisphere (excluding 
mountainous regions) (Pulliainen, 2010). The SWE prod- 
uct is based on the method of Pulliainen (2006). By incor- 
porating station observations and snow pack modeling into 
passive microwave retrieval algorithms, the goal is to pro- 
vide an accurate product useful for analyses at many differ- 
ent spatial scales, and for near-real time as well as 
climatological studies. The snow extent product is created 


A. Frei et al. I Advances in Space Research 50 (2012) 1007-1029 


1015 


using European Space Agency satellite visible and infrared 
observations (ERS-2 ATSR-2 and Envisat AATSR) based 
on the method of Metsamaki et al. (2005), and will likely be 
available at a variety of spatial resolutions. GlobSnow is 
currently available (http://www.globsnow.info/) but is 
new, so there is little peer-reviewed literature on it at the 
time of this writing (Takala et al., 2011). 

3.3.3. Other combined products 

A variety of combined products have been produced 
globally, regionally, or for specific purposes. One widely 
used combined product is NOAA’s National Operational 
Hydrologic Remote Sensing Center Snow Data Assimila- 
tion System, which operationally incorporates input from 
snow models, station reports, and airborne sensors to esti- 
mate daily SWE at 1 km resolution across the continental 
US (Carroll et al., 2001). The product by Brown et al. 
(2003), which employs the operational snow depth routine 
of the Canadian Meterological Center model (Brasnett, 
1999), has been used for evaluation of climate models (Frei 
et al., 2005). Foster et al. (2008) recently produced a global 
product blending visible and infrared, passive microwave, 
and active microwave scatterometer data, with the inten- 
tion of incorporating the most reliable aspects of each 
product. Derksen et al. (2004) produced a product going 
back to the early 20th century for the North American 
Prairies and Great Plains based on passive microwave 
and station observations. Biancamaria et al. (2011) esti- 
mated Northern Hemisphere fields of SWE based on pas- 
sive microwave combined with a dynamic snow grain 
model. Grundstein et al. (2002) developed a research-ori- 
ented SWE climatology for the Great Plains of the United 
States by combining station observations with the 1 -dimen- 
sional snow pack model SNTHERM (Jordan, 1991). A 
research-oriented product based on spatial interpolation 
of in situ depth measurements over North America (Dyer 
and Mote, 2006) has been used for process studies (Ge 
and Gong, 2008). The QuickSCAT active microwave scatt- 
erometer has been used to estimate the timing of snow melt 
across Greenland (Nghiem et al., 2001) and across Arctic 
lands (Wang et al., 2008). 

4. Methodology to compare and contrast products 

In this section we describe the methodology that we use 
to demonstrate the regions over which the products typi- 
cally differ. This analysis is not meant to provide insight 
into new remote sensing techniques, but rather to demon- 
strate the spatial extents and magnitudes of the differences 
between products during different seasons. The methodol- 
ogy employed here is designed to achieve two goals: (1) 
to identify regions across the Northern Hemisphere where 
there is agreement/disagreement between the three main 
products discussed here during clear days; and (2) to pro- 
vide an indication of the spatial distribution of uncertainty 
in the AMSR-E snow depth estimates, as determined by 
comparison to the CMC product. In this report we show 


results for three months: October (a month of rapid 
average increase of snow area), January (the month of larg- 
est average snow area), and April (a month of rapid aver- 
age decrease of snow area). For our analysis, AMSR-E 
SWE values are converted to depth. This is done using a 
fixed density mask, which is also used as part of the stan- 
dard product algorithm to estimate SWE values. We 
reverse the process in order to convert SWE values to 
depth. The reprojection methods, and the methods for each 
goal, are discussed in more detail in Sections 4. l^f.3. 

Many of our methods for goal 1 closely follow Frei and 
Lee (2010), and the reader is referred to that article for 
more details and justification of the methods. Note that, 
without independent verification, agreement between prod- 
ucts does not guarantee that they are correct; and, that if 
two of the three products agree, it does not guarantee that 
the third product is incorrect. 

4.1. Reprojection procedure 

IMS and MOD10C1 data sets were reprojected to the 
EAS E-Grid 25 km projection (Brodzik and Knowles, 
2002) (AMSR-E is already in this projection). Each 
EASE-Grid cell value was calculated as a binary (i.e. snow 
or no-snow) value. Because reprojection can introduce 
spurious errors at the grid-cell scale, and these errors 
are likely to be exacerbated in areas of variable terrain, 
we show no results for EASE-Grid cells within which the 
GTOPO 30 DEM elevation field has a standard devia- 
tion >100 m. We also avoid drawing conclusions from 
individual grid points, but rather focus on results over large 
regions with relatively little topographic variation. The 
reprojection, binary snow value calculation, and terrain 
masking were performed according to the method of Frei 
and Lee (2010). 

4.2. Agreement! disagreement between IMS, MODIS, and 
AMSR-E snow extent 

Since both IMS and MOD10C1 provide binary values 
indicating either the presence or absence of snow (the stan- 
dard MODIS products also provide fractional snow cover) 
but not snow depth, AMSR-E snow depths were converted 
to a binary value to facilitate this comparison. All AMSR- 
E depth estimates below 5 cm are considered snow free as 
that is the depth value assigned to shallow snow. 

All snow extent analyses include, at each grid cell, only 
days with “clear” skies, and only days for which all three 
products have non-missing data. We use the MOD10C1 
cloud mask to identify EASE-Grid cells that are mostly 
clear. Because MOD10C1 0.05 degree cells are higher reso- 
lution than the EASE-Grid and include fractional cloud 
cover, they can be used to estimate fractional cloud cover 
within each grid cell of our analysis. And, because the 
MOD10C1 cloud mask is considered conservative (in the 
sense that cloud-covered scenes are unlikely to be desig- 
nated as “clear”) (Riggs and Hall, 2002), we feel confident 
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that information from cloudy days is not being retained for 
analysis. This is achieved by retaining for analysis, for each 
EASE-Grid cell, only days with >80% of MOD10C1 cells 
that are <20% cloud covered, for which no product is miss- 
ing data. Frei and Lee (2010) present the rationale for this 
method and explain how the results are insensitive to rea- 
sonable values of these parameters. For each grid point 
on each day, either all three products agree (i.e. either snow 
or no-snow), or one product differs from the other two. 
The figures summarizing our results show, for each month, 
where each product disagrees with the other two products. 

4.3. Comparison of AMSR-E and Canadian Meteorological 
Centre (CMC) snow products 

For passive microwave data and the CMC model, no 
cloud mask is invoked, so we retain for analysis all available 
dates. While passive microwave data are not limited by 
most clouds, clouds with high liquid water content can 
affect the comparison between spaceborne- and ground- 
based SWE estimates (Wang and Tedesco, 2007); this issue 
is ignored here in order to increase the sample size. 

The comparison of AMSR-E to the CMC snow product 
is done by comparing climatological maps (2003-2010). 
For each month, three panels are shown containing maps 
of AMSR-E snow depth, CMC snow depth, and the differ- 
ence between the two products (we calculate the difference 
as CMC minus AMSR-E, so that a negative difference indi- 
cates that AMSR-E overestimates snow depth with respect 
to CMC). 

5. Results 

In this section we show the results of our analysis, the 
purpose of which is to demonstrate the spatial patterns 
of disagreement between the data sets. We also discuss pos- 
sible reasons for disagreements. In some cases these reasons 
may be speculative. 

5.1. Number of days per month available for analysis 

Before discussing disagreements between the products, 
we first show maps of the number of days per month avail- 
able for comparison (Fig. 4) which demonstrate the prob- 
lem presented by clouds. During October and January 
(Fig. 4a and b) most Arctic land surfaces are colored green 1 
or dark blue, which indicates that on average less than 
three (green) or three to six (dark blue) days per month 
are available for analysis. (In January one also sees the 
“ring” around the Arctic with no data associated with no 
solar illumination.) During spring, which tends to be less 
cloudy over most regions (Fig. 4c), one can find large por- 


1 For interpretation of color in Fig. 4, the reader is referred to the web 
version of this article. 


tions of the Arctic with either six to nine or nine to 12 days 
per month available for comparison. 

The vast majority of the unavailable days are caused by 
clouds, not by data that is missing for some other reason. 
Any satellite product based on visible and infrared band 
radiances will lack information from the surface under 
clouds. While passive microwave based products can pro- 
vide information under most types of clouds, they are 
currently unreliable under a number of circumstances (dis- 
cussed in the next section). Considering the importance of 
having daily real-time information about the surface to 
specify boundary conditions in weather prediction models, 
as well to track climatological changes in snow extent, 
IMS, or an equivalent product that provides information 
for all days regardless of cloud conditions, is a necessity. 

5.2. Disagreement between AMSR-E and the other two 
products 

Fig. 5 includes, for each month, a map showing where 
AMSR-E identifies snow to the exclusion of the both 
MOD10C1 and IMS (Fig. 5a, c, e) and a map showing 
where AMSR-E finds no snow when the other two prod- 
ucts identify snow (Fig. 5b, e, f). The most prominent fea- 
ture is the red colored plateau region of central Asia seen in 
all maps down the left hand column (Fig. 5a, c, e). This 
indicates that during all months over this region AMSR- 
E identifies snow more often than the other two products. 
While we do not, in general, assume that a product is 
wrong because it disagrees with the other two products, 
in this region we know from other studies that AMSR-E 
observations are biased due to problems in passive micro- 
wave snow detection at higher elevations associated with 
atmospheric influences on passive microwave retrievals 
(Wang and Tedesco, 2007). Since the atmosphere over 
the high elevation plateaus is much thinner, the algorithms 
calibrated globally at lower elevations require correction 
(Savoie et al., 2009). 

Panels on the right side of Fig. 5 show that during each 
month there are regions where AMSR-E identifies snow 
less frequently than MOD10C1 or IMS (Fig. 5b, d, f). 
The regions shown on these panels coincide with bound- 
ary of the continental snow cover during each month 
(see the Rutgers University Global Snow Lab web site 
for climatological maps of monthly snow cover based 
on IMS: http://climate.rutgers.edu/snowcover/index.php). 
Regions near the boundary tend to have patchy, shallow 
snow packs. During spring (Fig. 5f) the disagreement 
across well defined ablation bands at the southern bound- 
ary of the continental snow pack is also likely due to sig- 
nificant areas of melting snow with liquid water in the 
snowpack. 

5.3. Disagreement between IMS and the other two products 

Fig. 6 demonstrates that the most prominent circum- 
stance under which IMS disagrees with the other two 
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(a) October 


(b) January 
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Fig. 4. Average number of days per month at each grid point during which skies are clear and all three products are available (i.e. non-missing data) 
during (a) October; (b) January; (c) April. 


products is during the spring ablation period near the 
boundary of the continental snow pack (Fig. 6e). This 
result is in agreement with recent studies (Brown et al., 
2007, 2010; Frei and Lee, 2010) which find that over the 
last decade or so the timing of spring ablation over North 
America is later, by up to several weeks in the central 
Canadian Arctic, according to IMS in comparison to other 
observations. The reasons for these discrepancies, which 
are found during the entire spring ablation season (April, 
May, and June; May and June not shown here) over the 
boreal forest as well as the tundra, are not understood, 
but may be related to geographic factors such as the forest 
type and/or the presence of numerous lakes in this (Derk- 
sen et al., 2005a; Rees et al., 2006). Investigations into the 
cause of this problem continue. 


5.4. Disagreement between MOD10C1 and the other two 
products 

The most interesting example of MOD10C1 disagreeing 
with IMS and AMSR-E is found during autumn over 
Eurasia. During October over a broad, seemingly incoher- 
ent region of Eurasia, predominantly over Scandinavia and 
northern Europe, MOD10C1 often identifies snow when 
the other two products do not (Fig. 7a). However, this 
region is not as incoherent as it may seem, as it corresponds 
closely to the boreal evergreen needleleaf forest as defined 
by analysis of MODIS reflectance (Friedl et al., 2002). Dur- 
ing November (not shown) we find a similar pattern, except 
the differences are more extreme and concentrated more 
over Scandinavia. The eastern Eurasian region over which 
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(a) October AMSR-E snow (b) October AMSR-E no-snow 



(c) January AMSR-E snow (d) January AMSR-E no-snow 



(e) April AMSR-E snow 



A 

(f) April AMSR-E no-snow 
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Mountain 
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■ 30 - 40% 
20 - 30% 
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Fig. 5. Percentage of available days during which AMSR-E disagrees with IMS and MODIS. (a) October; AMSR-E shows snow while other two products 
show no snow, (b) October; AMSR-E shows no snow while other two products show snow, (c) January; AMSR-E shows snow while other two products 
show no snow, (d) January; AMSR-E shows no snow while other two products show snow, (e) April; AMSR-E shows snow while other two products 
show no snow, (f) April; AMSR-E shows no snow while other two products show snow. 


MOD10C1 often fails to identify snow when IMS and region of deciduous needleleaf forest. It seems that over 

AMSR-E see snow (Fig. 7b) corresponds closely to the one type of forest MODIS sees snow more often, while 
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(a) October IMS snow 


(b) October IMS no-snow 



(c) January IMS snow (d) January IMS no-snow 



(e) April IMS snow (f) April IMS no-snow 

Fig. 6. Same as Fig. 5 except showing where IMS disagrees with other two products. 


over a different type of forest MODIS sees snow less often. 
While the difficulties of remotely sensing snow under forest 
canopies have been widely reported, the patterns reported 
here have not been examined in the literature. 


5.5. Comparison of AMSR-E to the CMC snow product 

Maximum October snow depth values over the North- 
ern Hemisphere are ~20 30 cm (Fig. 8). The AMSR-E 
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(a) October MODIS snow 


(b) October MODIS no-snow 





(c) January MODIS snow 



(d) January MODIS no-snow 


(e) April MODIS snow 


(f) April MODIS no-snow 


Fig. 7. Same as Fig. 5 except showing where MODIS disagrees with other two products. 


product suggests more snow in Siberia than the CMC 
product. AMSR-E overestimation with respect to CMC 
over Siberia increases as the snow season progresses. In 
January, snow depth differences between the two products 


increase to ~30-40 cm (Fig. 9). In April, the area over 
which AMSR-E overestimates snow depth increases even 
further with respect to January. In contrast, over other 
regions AMSR-E tends to underestimate snow depth with 
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Fig. 8. Comparison of ASMR-E and CMC snow depth for October, (a) CMC estimated monthly mean snow depth, (b) AMSR-E estimated monthly 
mean snow depth, (c) Difference (CMC - AMSR-E) in monthly mean snow depths. 


respect to the CMC product, but these areas do not appear 
to expand as the snow season progresses. These include the 
Tibetan plateau and along the north-east coast of North 
America (Fig. 10). 

Histograms of the snow depth differences between the 
two products are shown in Fig. 11. Overall, AMSR-E tends 
to overestimate the values provided by CMC. While the 
variance of the errors can be seen in the histogram plots, 
perhaps a more informative number would be the coeffi- 
cient of variation (Cv). Cv defined as the absolute value 
of the standard deviation of the differences divided by the 
mean CMC snow depth, provides an indication of how 


large the differences are in comparison to the snow depth. 
For example, a value of Cv = 1 means that the errors are 
of the same magnitude as the mean depth; Cv = 0.1 means 
that the errors or an order of magnitude less than the mean 
snow depth. Cv values were calculated for each month 
(Table 1). Cv values are highest in October, when depths 
are small; lowest in January; and increase again in April. 
As a snow pack ages, even under cold conditions without 
additional precipitation, metamorphic processes lead to 
grain size variations (such as depth hoar formation) that 
tend to introduce errors in the passive microwave product. 
Furthermore, as temperatures fluctuate and additional 
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Fig. 9. Same as Fig. 8 but for January. 


precipitation events add fresh snow, snow packs can 
develop a series of well defined layers of different grain sizes 
that confound passive microwave based estimates of depth 
and SWE. Ice layers, which can develop as a result of 
melt-freeze and rain-freeze events, introduce additional 
scattering and therefore additional uncertainty. Such com- 
plications, combined with the impact of vegetation, 
especially vegetation that can change seasonally, can intro- 
duce a growing error in passive microwave retrievals as the 
snow season progresses. 

Improved confidence in our abilities to estimate snow 
mass from satellites would support efforts to monitor the 
fresh water flux into the Arctic Ocean. An order of magni- 
tude estimate suggests that the volume of water in the snow 


pack can play a significant role in the total annual river run- 
off into the Arctic Ocean of 4300M800 km 3 yr _1 (Arnell, 
2005; Miller and Russell, 2000). Our AMSR-E based 
(highly uncertain) estimate of the mean snow mass over 
land surfaces during March (the month of maximum snow 
mass) north of 60 N is ~1600 km 3 . Frei et al. (2005) based 
on the analysis of Brown et al. (2003) estimated the 
observed mean snow volume over North America during 
March to be ~1500 km 3 , which was equal to the median 
value estimated by a group of 18 climate models. This com- 
pares to a recent passive microwave-based estimate of 
~I400 km ’ and ~2300 km 3 for mean North American 
and Eurasian snow volumes, respectively (Biancamaria 
et al., 2011). The errors associated with most of these 
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estimates are currently unknown, but they indicate that the 
snowpack provides a significant fraction of the total river 
runoff to the Arctic. 

6. Discussion and conclusions 

For most of the snow season and most regions there is 
large-scale agreement amongst the products in identifying 
the location of snow covered surfaces (i.e. snow extent, 
regardless of snow depth) during clear sky conditions. 
One exception to this is over central Asia. It is known that 
passive microwave products identify snow on the Tibetan 
Plateau and surrounding mountains when visible and 
infrared products do not (Armstrong and Brodzik, 2001; 
Basist et al., 1996). Because passive microwave retrieval 


algorithms are calibrated at lower elevations, at these high 
elevations the reduced atmosphere between the surface and 
radiometer can result in retrieval errors (Savoie et al., 2009; 
Wang and Tedesco, 2007). The second exception occurs 
where snow is ephemeral, patchy, or wet. In such regions 
the attenuation of the passive microwave signal, upon 
which snow detection is based, is compensated for by emis- 
sion from the surface or from liquid water in the snow pack 
(Matzler, 1994; Ulaby and Stiles, 1980). Despite these dif- 
ficulties, all estimates (discussed in the preceding section) 
indicate that the snowpack is the source of a significant 
portion of runoff into the Arctic basin. 

The disagreements in snow extent during April are 
greater than during October or January in terms of the per- 
centage of available days during which one product differs 
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Fig. 11. Histograms of differences in snow depth (CMC - AMSR-E) for 
all grid points for (a) October, (b) January, and (c) April. 

from the other two. This is in agreement with Brown et al. 
(2010), Frei and Lee (2010), and Brown et al. (2007), who 
find differences between sensors during spring over North 
American regions experiencing ablation, and indicates that 
the wet snow during ablation is perhaps more of a hindrance 
to the identification of snow from satellites than some of the 
other confounding factors. However, during fall and winter 
the evaluation is hampered by data availability problems 
associated with cloudiness and solar illumination issues. 

Our analysis also demonstrates that snow depths esti- 
mated by the Canadian Meteorological Centre product 


Table 1 

Mean snow depth from CMC product; standard deviation of the 
differences between the CMC and AMSR-E snow depths; and the 
coefficient of variation. All values are averages of grid points across all 
Northern Hemisphere land areas north of 30 N excluding the Greenland 
ice sheet. 



(mean CMC 
snow depth) 
(cm) 

a (standard 
deviation of 
difference) (cm) 

Cv abs(CT/j) (coefficient 
of variation) ( unitless) 

October 

5.0 

4.15 

0.83 

January 

22.5 

7.18 

0.32 

April 

22.9 

12.33 

0.54 


and by the AMSR-E algorithm can differ substantially. 
Although there are no absolute surface reference observa- 
tions in most regions to determine which (if either) product 
is correct, we know from experience as well as theory that 
the passive microwave depth and SWE algorithms are inac- 
curate under certain conditions (Tedesco and Narvekar, 
2010). Sources of error include: surface heterogeneity 
within a passive microwave footprint; temporal and spatial 
variability in grain size and snow density; obscuration of 
snow by forests; masking of the passive microwave signal 
by liquid water in the snow pack; and effects of atmo- 
spheric attenuation. The persistent underestimation by 
AMSR-E with respect to CMC over some regions can be 
partially explained by considering that snow depth over 
many of those areas is above the ‘saturation’ depth to 
which the passive microwave algorithm is sensitive (Derk- 
sen, 2008; Markus et al., 2006; Matzler, 1994; Schanda 
et al., 1983); the presence of a high fraction of lakes over 
the north east of North America is also believed to be a 
source of error (Derksen et al., 2005a; Rees et al., 2006). 

Another example is the overestimation of snow depth by 
AMSR-E over northern Siberia, which can be attributed to 
the limitation of the current AMSR-E algorithm to 
account for the large grains that typically develop in snow 
packs in this (and some other) regions (Clifford, 2010). 
Over regions that develop and maintain a snow pack early 
in the season, the snow tends to insulate the ground keep- 
ing it warm even as air temperatures fall, resulting in a 
strong vertical temperature gradient in the snow pack. This 
temperature gradient causes vertical energy and vapor 
fluxes within the snow pack, the net effect of which is a 
layer of depth hoar at the bottom of the snow pack (Jordan 
et al., 2008). The large crystal sizes of depth hoar (~5 mm) 
cause increased scattering of microwave radiation resulting 
in an overestimation of the snow pack by the passive 
microwave algorithms. 

Opportunities remain for the development of improved 
snow products. For example, improvements can be made 
with regard to the retrieval of snow amount from passive 
microwave sensors (Tedesco et al., 2004) under forested 
terrain (Derksen, 2008), the refinement of snow extent esti- 
mates from visible and infrared sensors (Parajka et al., 
2008), and the estimation of sub-grid scale information. 
Tedesco and Miller (2007b) explore the relative merits of 
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combining active and passive microwave retrievals, using a 
MODIS snow product as their reference “truth.” A num- 
ber of researchers are investigating the potential for finer 
scale information on snow extent, amount, fractional snow 
cover (Derksen et al., 2005b; Salomonson and Appel, 2004, 
2006), snow melt (Wang et al., 2008), as well as on snow 
pack properties (Kinar and Pomeroy, 2007; Nolin and 
Dozier, 2000; Painter and Dozier, 2004; Painter et al., 
2003; Rango et al., 2000; Schmugge et al., 2002). Improve- 
ments in remotely sensed products that do not rely on the 
assimilation of data or model results will come as a conse- 
quence of improved understanding of the interaction 
between electromagnetic and geophysical parameters at 
large spatial scales. In this context, a new operational algo- 
rithm based on the inversion of an electromagnetic model, 
artificial neural networks and snow climatology currently 
under evaluation may be capable of accounting for some 
of these limitations. 

One currently active area of research is the development 
of combined products, which include in situ observations 
and/or modeling results as well as remotely sensed infor- 
mation. One can identify advantages and disadvantages 
to both combined and stand-alone remotely sensed prod- 
ucts. While stand-alone remotely sensed products contain 
inherent drawbacks as discussed here, at any time, either 
in situ or remotely sensed data streams can fail, rendering 
combined as well as stand-alone products vulnerable to 
missing information. This is most critical for real- or 
near-real time operational products, on which weather 
forecast models or time-sensitive decisions rely. 

Remote sensing of snow continues to contribute to our 
understanding of Earth system processes. MODIS snow 
products are valuable because they can provide high reso- 
lution snow estimates under cloud-free conditions using a 
quantifiable algorithm. However, for climatological as well 
as operational purposes, humans can integrate and filter 
data from multiple sources and satellite images in ways 
that fully automated methods are (at least currently) 
unable to, and provide information for the entire land sur- 
face of the globe, regardless of the presence of clouds. 
Thus, continuation of IMS, with its long record of snow 
extent, is a priority. Considering the difficulties in deter- 
mining SWE on a global scale from stand-alone remote 
sensing products, it seems likely that combining multiple 
sensors with station observations and/or models, such as 
in the GlobSnow product will provide the best estimates 
of SWE. 
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